Coherent electron-phonon states in suspended quantum dots: decoherence and 

dissipation effects 
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The dynamics of coherent electron-phonon (el-ph) states is investigated for a suspended nanos- 
tructure. Exact quantum dynamics calculations reveal that electron and phonons (comprising a 
thermal bath) couple quantum mechanically to perform coherent oscillations with periods in the 
range of tens of nanoseconds, despite the finite temperature of the phonon bath. Mechanical energy 
dissipation due to clamping loss is taken into account in the calculations. Although the lifetime of 
the coupled el-ph states decreases with the temperature, well defined Rabi oscillations are obtained 
for temperatures up to 100 mK. The dynamics of the coupled electron-phonon state is susceptible 
to various forms of external control. For instance, a weak external magnetic field can be used to 
control the dynamics of the system, by decoupling the electron from the phonon bath. The results 
cast light upon the underlying physics of a yet unexplored system that could be suitable for novel 
quantum device applications. 
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INTRODUCTION 

The possibility of engineering semiconductor devices 
at the nano and micro scales has created the conditions 
for testing fundamental aspects of quantum theory oth- 
erwise difficult to probe in natural atomic size systems. 
Particularly, quantum dot (QD) systems have been rec- 
ognized as a physical realization of artificial atoms and 
molecules, whose properties (e.g., structural and trans- 
port) have been intensively investigated, notably in the 
presence of magnetic fields [H, HI & HI HI • QDs have be- 
come the building blocks of various quantum devices and, 
nowadays, they can be coupled in arrays to create charge 
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gates. 



A new promising possibility for both implementing and 
investigating coherent phenomena in semiconductor de- 
vices is the combination of q uantum dot systems with 
suspended nanostructures [ll[ 12, 



13, 14, 15, Id. As an 



immediate consequence, it is possible to improve the iso- 
lation of the electronic quantum system from the bulk of 
the sample. But it is also expected that nanoelectrome- 
chanical systems (NEMS) will lead to the investi gati on of 



new regimes of phonon-mediated processes [111 1 1 61 ] and 
the observation of quantum behavior in mesoscopic me- 
chanical systems. In addition, phonon cavities can be 
envisaged as a solid state analog of quantum electrody- 
namic cavities. 

For all these pursuits it is generally desirable to con- 
struct systems with very little loss of energy and very 
high quality factors Q. High frequency resonators have 
been fabricated from a variety of materials, most of them 
based on silicon. Due to potential applications as on-chip 
high-Q filters and clocks, nanomechanical resonators are 
rapidly being pushed to smaller size scales and higher 
frequencies. Resonant frequencies higher than 1 GHz 
have consistently been achieved 17|, [l8(. However, it 



has also been observed that the quality factor Q of mi- 
cro and nano resonators decrease as the size of the sys- 
tem is reduced and the fundamental frequency increased 
[H HI H3. That effect imposes severe constraints on 
the design of NEMS oscillators, specially for those ap- 
plications that rely on the response of a single resonant 
frequency. Therefore, it is important to look for alterna- 
tive systems that can preserve coherence even in the case 
of high mechanical dissipation. 
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FIG. 1: Illustration of the NEMS structure comprised of a 
circular quantum dot of radius R — 75 nm embedded on a 
suspended platform with dimensions w — 1.2 /im (width), 
I — 200 nm (length) and t — 50 nm (thickness). The figure 
also depicts a metallic gate extending through the bulk of 
the sample, which affords external control of the state of the 
electron with minimal interference with the vibrational state 
of the platform. 



Here we investigate the underlying physics of a novel 
type of nanoelectromechanical oscillator that is based on 
the dynamics of coupled electron-phonon states. The 
work looks into the evolution of a QD state with well de- 
fined electronic angular momentum as it interacts with 
a thermal field comprised by the phonon states of the 
platform at finite temperature. As a consequence of the 
electron-phonon interaction, the electronic angular mo- 
mentum oscillates between degenerate eigenstates of L Zl 
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bridge structures. The most serious difficulty is the 
augmented dissipation of vibrational energy through the 
wide clamping edge of the platform. However, we find 
that the deleterious effects caused by the low Q factor of 
the platform are compensated by the fact the dynamics 
of the QD is determined from outset by its interaction 
with a thermal field. The influence of external and inter- 
nal mechanisms of energy dissipation will be discussed in 
the light of the simulation results. 

Section presents the theoretical methods developed to 
calculate the quantum dynamics of the proposed nano- 
electromechanical structure at finite temperatures. Sec- 
tion is dedicated to the analysis of the simulation results 
under different physical conditions. The conclusions and 
a final discussion are presented in section . 



THEORETICAL MODEL 



FIG. 2: Contour plots of the transverse vibrations of the plat- 
form. From the bottom to the top, the first 5 phonon modes 
are depicted. The phonon modes have even (+) or odd (- 
) parity with respect to the width dimension. The letters C 
and F denote the clamped and free edges. The circle represent 
the position of the QD. 



like in a resonant Rabi system. A weak magnetic field 
can be used to decouple the electron from the thermal 
bath, by splitting the L z eigenstates without disturbing 
the original symmetry of the electronic wavefunctions. 

For that purpose, we consider a yet unexplored design 
for the suspended NEMS, whose width (w) is larger than 
the length (7), as shown in Figure [TJ For the sake of con- 
venience we name it platform, to make the distinction 
with the usual cantilever structures. The dimensions of 
the suspended platform used in the following calculations 
are w = 1.2 fim, I = 200 nm and t = 50 nm (thick- 
ness). Among the reasons that moved us to consider 
such NEMS, we point out that its fundamental frequency 
(/i ~ 2 GHz) is approximately 40 times higher than 
the fundamental frequency of a cantilever with equivalent 
dimensions (i.e., w = 200 nm, I — 1.2 fim). The short 
spaced density of states of a normal cantilever (or bridge) 
favor dissipation and wash out the coherent electron- 
phonon dynamics. In addition, the proposed configu- 
ration allows a direct external contact with the electron 
in the quantum dot, by means of an electrostatic gate 
extending through the bulk of the sample, as depicted in 
the Figure Q] We notice that such a noninvasive exter- 
nal coupling with the QD is not possible in the regular 
cantilever or bridge structures. The lowest transverse vi- 
brational modes of the platform are depicted in Figure [2] 
together with the circular QD. 

Despite the advantages, the proposed NEMS poses sig- 
nificant caveats compared to the usual cantilever and 



Suspended electron-phonon system 



The hamiltonian of the system is written as 
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— H e i + H p h + H e i_ p h + H re i ax 
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(1) 



which includes the electron and phonon hamiltonians 
(Hq = H e i + Hph), the electron-phonon interaction 
(H e i_ p h) and the phonon relaxation term (H re i ax ). In 
brief, b\ and b K are the electron creation and annihila- 
tion operators, while a a and a a are the corresponding 
phonon operators, with h a = a^a Q . The el-ph matrix el- 
ements are denoted by g K > aK . The energy relaxation rate 
of the phonon system with the external environment is 
7 Q = huj a /Q, where Q is the quality factor of the plat- 
form, rTa~(T) is the thermal occupation of phonon mode 
a and (n a ) the corresponding quantum average for the 
coupled el-ph state. 

Hamiltonian (jTJ) is a generalization of the usual atom- 
field interaction [22j |. here describing a multi-level elec- 
tronic system in a multi-mode bosonic cavity . In the 
following calculations we assume that the QD is occupied 
by a single electron - such devices have been extensively 
characterized in the literature @, [H, [H [H, [26j|. The 
model can also describe the excess electron of a QD. The 
one electron hamiltonian for a quasi two-dimensional cir- 
cular QD of radius R has eigenstates 



<Pk(T,9) = 
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with k = (I, v), I — 0, ±1, ±2, . . . and ai v the v-lh root 
of the Bessel function of order |Z|, Jm (pn v x). The cor- 
responding energies are E K — ^-j^t, where m e is the 
effective mass of the electron. An interesting condition 
arises if a weak magnetic field B is applied perpendicu- 
larly to the plane of the quantum dot 2jj • In this case 
the one electron hamiltonian becomes 
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where V(r) is zero for r < R and infinity outside the 
dot. The vector potential A = B/2 (—y,x,0) is written 
in the symmetric gauge. For weak magnetic fields, such 

1 /2 

that the magnetic length Is = (hc/eB) 1 is larger than 
the radius of the dot, R < Ib, the diamagnetic term 
{Ediam = $ m ^ c i r 2 ) in Eq (g]) can be disregarded and the 
one electron energies can be written as 



Eiu — 
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+ fi B mB 



(5) 



where fi B — eh/2m e c is the Bohr magneton. Moreover, in 
the weak field approximation the wavefunctions ([2]) are 
still eigenstates of H e i, since L z — —ihd/dO commutes 
with the hamiltonian. For R = 75 nm and B = 500 G, 
such that l B — 1.5R, E diam / Ezeeman = 0.1 and the error 
committed by Eq. ([5]) in the calculation of the energy is 
less than 1% (more details in Section ). 

The phonon eigenstates of the platform are obtained 
through the quantization of its mechanical modes of vi- 
bration, according to the formalism detailed in Ref [ll| . 
The operators and a a , associated to h a = aJ,a Q in Eq. 
(fTj), create and annihilate phonon modes in the cavity. 
At low temperatures only the long wavelength acoustic 
modes are relevant and the semiconductor is described 
by the elastic approximation, so that the classical plate 
theory (CPT) [28| can be used to calculate the vibra- 
tional modes of a thin platform. The fundamental fre- 
quency calculated for the suspended platform depicted 
in Fig. [T] is /i « 1.88 GHz. The dynamics of the cou- 
pled electron-phonon bath state depends on the spectrum 
of the phonon modes, therefore it is important to vali- 
date the calculated phonon energies. For that purpose 
we have compared the CPT modes with those obtained 
from the rigorous three-dimensional (3D) elastic equa- 
tions [2!|, [301 • It is found that the CPT method consis- 
tently yields the correct mode shapes and, moreover, the 
fundamental frequency /i evaluated by the CPT method 
is just 3% higher than the value yielded by 3D method. 
It is also well known that the CPT method tends to over- 
estimate the vibrational frequencies as the order of the 
modes increase, however, the energy discrepancy of the 
higher modes is compensated by their low thermal occu- 
pation for temperatures T < 100 mK. 



The hamiltonian for the electron-phonon interaction 
reads 



H, 
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where g K >aK, = CdpV®^^/ y/U^, in connection to Eq. 
(U), Cpd denotes the deformation potential (DP) con- 
stant of the material and b\. (b K ) is the fermionic cre- 
ation (annihilation) operator satisfying the usual anti- 
commutation relations. The electron-phonon interaction 
is written in terms of the matrix elements V K ' QK that de- 
scribe the coupling between electrons and phonons in the 
NEMS. The properties of V K / aK depend on both the mate- 
rial attributes and the geometrical symmetries of the de- 
vice. At low temperatures the more important electron- 
phonon coupling mechanisms are described by the defor- 
mation (DP) and the piezoelectric (PZ) potentials, how- 
ever, only the former is included in ([1]), because silicon is 
not a piezoelectric material. A detailed derivation of Eq. 
([6]), including the piezoelectric interaction, can be found 
in Ref [16(. Furthermore, it has also been demonstrated 
[ill that, for a circular QD located at the center of the 
platform, the even parity phonon modes produce real val- 
ued matrix elements g K ' aK whereas the odd parity phonon 
modes yield imaginary valued elements. If the QD is dis- 
placed from the center of the platform g K i aK is complex. 
Moreover, the matrix representation of the total hamil- 
tonian H can be written in blocks of fixed electronic an- 
gular momenta Vs as H;'; = Ain + (i) 



sign(i'-Z) 



H'l 



where 



A and B are real, symmetric and mutually disjoint matri- 
ces. Taken individually, the matrix A (B) is responsible 
for the coupling between the electron and the even (odd) 
parity phonon modes. 

Finally, the dominant source of mechanical energy dis- 
sipation in the NEMS is due to attachment loss via the 
platform clamping edge. The mechanism is described 
within the relaxation time approximation in the hamil- 
tonian (fT]), with 7 Q = fkv a /Q, where Q is the quality 
factor of the NEMS, assumed to be the same for all 
the modes. In Eq. [j] n^(T) = l/(exp[fuv a /k B T} - 1) 
is the thermal occupation of the phonon mode a and 
(n a ) — Tr{n a p(t)} denotes the quantum average of the 
phonon occupation operator. The later depends on the 
time through the electron-phonon density matrix p(t). 
Then, if (n Q ) > n^(T) during the time evolution of the 
system, because of the energy that is exchanged with the 
electron, mechanical energy flows out of the plataform to 
the bulk of the sample. The relevance of other dissipation 
mechanisms is discussed in section . 

The electron-phonon basis set is formed by the di- 
rect product of the one-electron states \4> K ) with the 
multi-phonon states |n) = |ni, ni, ns, tin). Here, 
n a = 0,1, ... ,n denotes the number of phonon quanta 
per mode a. A total of N distinct phonon modes are 
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considered, so that a typical basis vector is written as 



N 



c;n) = \<j> K ) <g> \ ) 



(air- io) 



(7) 
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with TV = 40 and n a ^ 40 in the following calculations. 
The convergence of the results were tested with respect to 
both the electron and phonon basis sizes. The electron- 
phonon basis is set up according to the prescription: first, 
a large (> 1.5 x 10 5 ) basis set comprised by states of Eq. 
([7]) is generated and energy sorted, in order to create 
a micro-canonical ensemble. The time evolution of the 
initial el-ph state is then carried out on a truncated basis 
with size ranging from 1 x 10 5 (for calculations made at 
T = 50 mK) to 1.2 x 10 5 (for calculations made at T = 
100 mK). For the parameters used in the calculations, 
the el-ph basis is formed by the combination of 50 to 60 
distinct electronic states \</> K ) and approximately 15 x 10 3 
different Ji2, na, ■■•,tin) phonon states. We point out 
that the calculations do not converge if the electronic 
basis does not include several angular momentum states. 
We have used 7 < |/| < 10, but the appropriate number 
depends on the characteristic energies of the system. 



Dynamics by the split-time Chebyshev method 

The time evolution of the coupled electron-phonon 
state is governed by the time-dependent Schrodinger 
equation 



ih 



dm 

dt 



H\9) , 



(8) 



with H given by the hamiltonian in Eq. ([!]). Within a 
sufficiently small time slice the time dependence of the 
hamiltonian can be disregarded and the solution of Eq. 
(JSJ) for a small time step At is 



|#(Af)) = exp 



|*(0)> = U(M) |*(0)) . (9) 



The time propagation of has been checked for de- 

creasing time steps in order to validate the approxima- 
tion. Within a given time slice, a very efficient propaga- 
tion scheme consists of expanding the evolution operator 
U in terms of the orthogonal Chebyshev polynomials T k 

IE El 
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with the coefficients c& (i) calculated using the orthogonal 
property of the Chebyshev polynomials 



Ck(t) = / — dx 

K J-l \/l~x 2 

= (2-S k )(-i) k J k (t) , 



(11) 



where Jfe(i) is the Bessel function of order k. Because 
the Chebyshev polynomials T/~ (x) = cos [k arccos(x)] are 
defined in the interval [—1,1], the hamiltonian H has to 
be rescaled according to 



H = 2 



(H-E) 



(12) 



where E = (E max + E m i n )/2 corresponds to the medium 
value of the eigenvalues of the original hamiltonian and 
e ^ E m ax — E m i n is a superior value for the range of its 
span, so that —1 $C (k; n\H \n; n) ^ 1. 

Following transformation (|12|) . the time step is rede- 
fined as t — (eAt/2h). Therefore, the time evolution for 
the initial state \^>(to)) = \^>o) is given by 



|*(At)) = e~ lHT e - lEAt / h \^ c 



(13) 
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with Cfe(r) = (2 — <5fc)(— i) fc Jfc(r). The vector states 
= Tfe(W)l^'o) ar e easily obtained from the recur- 
rence relation 



2H|* fc _!) - |$ fc _ 2 



(15) 



and from the particular results Tq(H)\ 1 ^o) = \^o) and 

T 1 (W)|* )=W|*o}- 

The advantage of the Chebyshev propagation method, 
compared to other expansion schemes is due to the fact 
that when k > r, Jfc(r) tends to zero exponentially fast. 
Therefore, expansion (TT4| can be truncated at K > t. 
with the actual value of K depending on the chosen ac- 
curacy e. In our calculations e — 5 x 10~ 5 , which yields 
K/t w 1.1 for a time step At = 0.25 ns. In the Cheby- 
shev method, CPU time scales linearly with basis size 
and memory use is highly optimized because the matrix 
elements of Ti. can be calculated on the fly, or stored in 
packed format for sparse hamiltonians. 



SIMULATION RESULTS 

The quantum-mechanical equations of motion of a two- 
level atom coupled to a single resonant mode of a cavity in 
the absence of dissipation were solved exactly by Jaynes 
and Cummings [331 ] . They found that this interaction 
leads to a simple oscillation of energy between the atom 
and the cavity mode. This result corresponds to the well 
known semiclassical solution of a two-level system inter- 



acting with the electromagnetic field (22|. Sachdev [34 1 
examined the quantum theory of a two-level atom in a 
damped cavity, in the limit of small damping, and showed 
that thermal fluctuations wash out the Rabi oscillations 
as the temperature increases. 

However, we find that a simple two-level model for 
the electron in the QD is not adequate to describe the 
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FIG. 3: Diagram with the energies involved in the problem: 
the energies Ei v of the eigenstates of H e i (see Eq.([5}) are la- 
beled by the angular momentum (I) and the radial quantum 
number (v); the energies of the lowest 40 phonon modes, and 
the thermal energy parameter (fc_gT). For the case of a per- 
pendicular magnetic field, the Zeeman splitting is indicated 
for the electronic states with 1^0. The magnetic field is 
B = 500 G. 



present results, even if the higher energy angular momen- 
tum states are not significantly populated. The reason 
is that the angular momentum of the electronic system 
is not stationary anymore, because of its coupling with 
the vibrational modes of the platform. Thus, although 
the initial state of the electronic system may consist of 
a wavepacket of well defined angular momentum, several 
other angular momentum states will take part in the dy- 
namics. 

In the remainder of the paper we investigate the co- 
herent dynamics of the electronic state of the QD cou- 
pled to a thermal field comprised by the phonon modes 
of the platform at the initial temperature T . The to- 
tal vibrational energy of the multi-mode cavity state 
n = (ni,n 2 , ■ ■ ■ ,n N ) is E ph (n) = J2 a (n a + |) tiUa, thus 
the density matrix elements of the phonon ensemble in a 
thermal field configuration are given by 



exp{-E ph (n)/k B T) 
£ {n} exp(-E ph (n)/k B T) 



(16) 



In the following calculations two cases are considered: 
the initial temperature of the thermal field is either T = 
50 mK or T = 100 mK, which correspond to typical 
conditions in the study of quantum coherent phenomena 

ansa. 

The relevant energy scales of the system are shown 
in Figure El where we have the phonon energies E p h ~ 
8 x 10~ 3 — 10 _1 meV, the electronic energies E e i ~ 
10 _1 meV, the electron-phonon matrix element energies 
Eei- P h ~ 10~ 7 — 10~ 4 meV and the thermal energy 
Eth ~ 10~ 3 — 10~ 2 meV. The Zeeman splitting between 
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states with I = ±1 is also shown, AEzeeman ~ 6 x 10 
meV for B = 500 G. The quality factor of the platform 
is Q = 100 for all phonon modes. The calculations are 



performed for a suspended platform made of crystalline 
silicon [35], with dimensions w = 1.2 fim, I = 200 nm 
and t = 50 nm and a circular QD of radius R — 75 nm. 

The el-ph dynamics is sensitive to the position of the 
quantum dot in the phonon cavity, due to the interplay 
between the distinct symmetries of the circular QD states 
and the rectangular phonon modes of the platform. In 
the following calculations, the QD is placed at the center 
of the platform, which is the most symmetric position. 



Intrinsic electron-phonon bath dynamics 

The dissipation mechanisms can be generally under- 
stood as having external or internal origin. Initially, 
we investigate the physics of the intrinsic el-ph bath re- 
laxation dynamics in the absence of magnetic fields for 
the idealized case Q — > oo. Because the electronic an- 
gular momentum L e i is stationary for the isolated QD, 
its time evolution is particularly useful in describing the 
electron-phonon interaction. Thus, we look into the time 
evolution of the rotational wavepacket \ty(t — 0)) — 
1 1,1) ® \xbath), which describes the single electron QD 
in the first excited state (L e i = 1) in contact with the 
phonon bath described by Eq. (flTJl) . 
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FIG. 4: Time evolution of the coupled electron-phonon state 
assuming only internal decoherence and no external mechan- 
ical dissipation. From top to bottom the graphs corre- 
spond to: A) average electronic angular momentum L e i(t) = 
Tr e i{Lp e i(t)}, B) electronic decoherence parameter T e i(t) = 
Tr{/5g Z (t)} and C) electronic energy. The curves correspond 
to T = 50 mK (black) and T = 100 mK (red). 

Allowing the initial state |W(0)) to evolve according 
to Eq. ([Tip. Figure HJ-A shows the electronic angu- 
lar momentum L e i(t) = Tr{p e i(t)L} as it oscillates in 
time between the degenerate states L e i — ±1, due to 
the interaction with the phonon bath in the platform. 
Pei{t) = Tr p h{p(t)} is the reduced electronic density ma- 
trix and Tr p h designates the trace over the phonon states. 
The black curves represent the calculations made for the 
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initial phonon temperature T = 50 mK and the red 
curves correspond to T = 100 mK. 

Rabi's theory is useful in interpreting the general fea- 
tures of the numerical results and can be used to de- 
termine the effective electron-phonon coupling parame- 
ter 7 ep responsible for the L e i oscillations. Analysing 
the dynamics of the simpler case, i.e. T — 50 mK 
in Fig. EJ-A, we observe that L e i oscillates with al- 
most constant amplitude between the values -1 and +1. 
In this case, L e i is equivalent to the inversion ampli- 
tude W(t) = E„(l^+i,n| 2 - |C-i,n| 2 ), where |tf(t)) = 
J2 K ,n C K,n\K,n}. Moreover, for T = 50 mK, (n a ) < 1 
for all phonon modes and one can assume the vacuum 
state for the phonon cavity. Therefore, L e \ can be ap- 
proximately described by 



J ±i 



L e l(t) = 



cos 



'±1 



<4l + (27e P A) 5 



,(17) 



where ui±i = (-Ei l — E_xi)/h. In the absence of mag- 
netic fields lj±x = and Eq. (jTTJ) simplifies to L e j = 
cos(2j ep t/h). The evaluation of the effective electron- 
phonon coupling yields j ep « 5.10 -5 meV for this sim- 
ple case, which is consistent with the calculated ma- 
trix elements <? K ' aK in Eq. [1] It is important to notice 
that the dynamics of the el-ph system is caused mainly 
by virtual phonon transitions. When a weak magnetic 
field is applied, the electron-phonon matrix elements do 
not change, but the degeneracy condition is broken and 
uj±i = 2Ezeeman/f l - Then, one expects that the oscil- 
lation amplitude of L e i decreases. Real transitions can 
also occur as a result of accidental resonances involving 
the / = state and high energy phonon modes. Those 
effects are discussed in section along with the simulation 
results. 

A variety of different behaviors is observed for L e i, just 
by changing the edge conditions of the platform or by dis- 
placing the QD from its center. The effects include the 
plain oscillatory dynamics, beatings and the overdamped 
decay of L e i . The asymptotic value of the electronic an- 
gular momentum for all the cases is L e i = 0. That cor- 
responds to the equal occupation of the electronic states 
with I = ±1, a small occupation of states with I = and 
a very small occupation of the other I states. 

As expected, the damping of the coupled electron- 
phonon oscillations varies with the initial bath tempera- 
ture. To evaluate more precisely the degree of decoher- 
ence in the electronic wavefunction (i.e., the rotational 
wavepacket) we calculate the purity of the electronic state 
T e i(t) = Tr{/5g Z (t)}. At time t — the electronic system 
is in the pure state described by |\&(t = 0)) and T e ; = 1. 
The electronic decoherence is presented in Figure 0J-B, 
for the initial phonon temperatures T = 50 mK (black) 
and T — 100 mK (red). The interaction with the phonon 
modes of the platform is responsible for the decay of T e i , 



which saturates at r e ; w 1/2, evidencing the fact that 
the phases of the wavepacket evolve towards an even dis- 
tribution among the \k; n) states with angular momenta 
I = ±1. The calculations show that the internal el-ph 
friction causes electronic decoherence in less than 100 ns 
for T = 100 mK. We point out, however, that such time 
scale for the electronic decoherence should be understood 
as a higher bound , since the external decohering mech- 
anisms, like the clamping loss, have not yet been taken 
into account. 

By examining the autocorrelation of the entire 
electron-phonon bath state, £(i) = (^(0)1^(4)), we have 
evidenced for all the cases investigated that the electron- 
bath system undergoes a fast phase decoherence, within 
the time scale T e i- P h ~ 0.1 ns. It is followed by the elec- 
tronic decoherence, evinced by the decay of r e ; within 
T e i ~ 100 ns i an d by an even slower internal energy dissi- 
pation process (el-ph friction) that leads to the formation 
of the polaron state in the QD. The electronic energy 
dissipation is shown in Fig0]-C, a stationary energy dis- 
tribution is reached within the microsecond time scale 
(r erg ). The initial energy of the electron, indicated by 
the arrow, is much higher than the thermal energy of 
the phonon bath, leading to a recurrent energy exchange 
between the electron and the phonon modes of the plat- 
form. Energy is exchanged between the electron and the 
phonon bath primarily through the I = state and the 
high energy phonon modes (a = 15, 16, 17). For the sake 
of comparison, similar beating effects have been observed 
in molecular dynamics simulations of nanomechanical en- 
ergy exchange between single-walled carbon nanotubes 



Dissipation effects 

Energy dissipation is a central issue for the operation of 
nano- and microelectromechanical structures (MEMS). 
The identification of the source of mechanical energy dis- 
sipation is, nonetheless, a complex problem because dif- 
ferent mechanisms can be determining in distinct phys- 
ical situations 0, [43[. Usually, the investigation is fo- 
cused on the relaxation of energy from a driven resonant 
mode of a cantilever, bridge or structure alike. For a sus- 
pended NEMS operating in vacuum, vibrational energy 
is radiated out of the structure through the resonator 
attachments 



21, 37, 38. 3 



Oj. On the other hand, 
intrinsic dissipation also contributes to take energy out 
of the resonant mode. Such mechanisms include phonon- 
phonon and electron-phonon scatterin g, th ermoelastic ef- 
fect [Ilj], surface and bulk defects [Hl2Q|, SI] . 

At sub-Kelvin temperatures and in the mesoscopic 
regime, phonon-phonon scattering and thermoelastic re- 
laxation can be safely disregarded 4l| in comparison to 
the other intrinsic friction mechanisms. Dissipation due 
to defects, however, can not in principle be discarded. 
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Defects arise most frequently from broken and dangling 
bonds on the surface of the structure and from contami- 
nation by other atoms. Zolfagharkhani and collaborators 
(i^ l have conducted a detailed investigation of the quan- 
tum friction in nanomechanical oscillators at millikelvin 
temperatures. They showed, in accordance with other 
studies (20j . that dissipation in this regime is dominated 
by the interaction of the resonant mechanical modes with 
localized defects. They also found that dissipation de- 
creases with temperature, but saturates at<3 _1 ~3T0 _5 
for T < 100 mK, in the single-crystal Si nanobridges in- 
vestigated. 

Undoubtedly, the dominant dissipation mechanism in 
the platform is attachment loss, because of its wide 
clamping edge. Throughout the literature, experimental 
and theoretical investigations concentrate on cantilever 
and bridge geometries, owing to the fact that the at- 
tachment losses can be minimized for high aspect ratio 
structures ( i.e., I 3> w) 21, 4(3]. Therefore, it is difficult 
to estimate the quality factor of the suspended platform. 
However, two different theoretical expressions yield the 
same estimate for the Q of our model structure. Ac- 
cording to the expression of Photiadis and Judge [3§| 
Q P j ps 3.2 (l 5 /wt 4 ) 137 whereas the estimate due to 
Jimbo and Itao ^ yields Qji « 2.17 (l/tf ps 139. Both 
estimates are in qualitative agreement with the experi- 
mental results of Huang and collaborators [2l| , obtained 
for a series of SiC nanobridges. That work measured a 
quality factor Q = 500, for a wide (w = 120 nm) and 
short (1 = 1 [im) doubly clamped beam, with funda- 
mental resonant frequency fi ~ 1 GHz. 

Therefore, we estimate that the overall quality factor 
of the platform of Fig. [I] is Q « 100. Assuming that 
the structure has Q=100, one finds that the lifetime of a 
resonant mechanical mode with frequency / = 1.5 GHz is 
limited by the losses in the platform to r Q « Q/2irf « 10 
ns. However, the dynamics of the el-ph system is caused 
by virtual phonon transitions that do not result in direct 
energy transfer to the phonon system, nor to a particular 
resonant phonon mode. In that case, one expects that the 
electronic oscillations should persist even in the presence 
of mechanical dissipation. 

Since the phonon system is described by a thermal field 
from the outset, the dissipation of the mechanical energy 
of the platform is described within the relaxation time 
approximation in hamiltonian ([1]), with Q — 100 for all 
phonon modes. According to H re i ax , mechanical energy 
of a given phonon mode is radiated out to the bulk of the 
sample if (n a ) > n^(T), and vice- versa. As evinced by 
Fig. [5l the mechanical dissipation produces a strong ef- 
fect on the low temperature el-ph bath dynamics (T = 50 
mK) whereas the T = 100 mK case is just weakly af- 
fected. The reason is that a small amount of energy that 
is transferred to the low temperature phonon system ren- 
ders enough to drive that system out of the equilibrium. 
The result also indicates that electron-phonon friction 




o. 



FIG. 5: Time evolution of the coupled electron-phonon 
dynamics, including attachment loss dissipation with Q = 
100. A) average electronic angular momentum L e i(t) = 
T<Cei{Lp e i{t)}, B) electronic purity T e i{t) = Tr{p 2 el {t)}. The 
curves correspond to T = 50 mK (black) and T = 100 mK 
(red). 



becomes the main source of electronic decoherence for 
T > 100 mK, but below that temperature clamping loss 
is determining. 



Dynamics of L e ; in the Bloch sphere 

The entanglement between the orbital angular momen- 
tum states |Z = ±l,i/) is revealed through the anal- 
ysis of the el-ph hamiltonian matrix i7 e /_ p /j. It has 
been noted that, for a circular QD located at the center 
of the platform, the matrix representation of the total 
hamiltonian H can be written in blocks of fixed I's as 



H'l 



sign(/'-/) - 



where A and B are real, 



symmetric and mutually disjoint matrices. Moreover, 
the matrix A (B) is responsible for the coupling between 
the electron and the even (odd) parity phonon modes. 
Therefore, starting from |*(i = 0)) = |1,1) <g> \xbath), 
the A matrix operates like the generator of rotations R x , 
leading the electronic state through a precessional move- 
ment around the e x axis of the angular momentum space. 
Likewise, the B matrix functions like the R y operator, ro- 
tating the electron state around the e y axis. 

Since the time evolution of |W(0)) results mainly in oc- 
cupation of the I = ±1 states with v = 1 we can, for 
the sake of clarity, restrict ourselves to that subspace 
and describe the electronic state in the Bloch sphere. 
An the arbitrary state in the Bloch sphere is written as 
= cos 1 1 + 1) 2 +exp(£<£) sin 1 1 - l) z , with < 6 < tt 
and < (j> < 2ir. In the current model the dynamics of 
|\&(f)) in the Bloch sphere is evidenced by the density 
matrix element p e ;(+l,— 1) = /9* ; (— 1, +1), with v = 1, 
whose modulus (\p e i(+l, — 1)|) and phase ((f)) are shown 
in Figure [5] for T — 100 mK and Q=100. The numerical 
calculations show that the electron state undergoes a pre- 
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FIG. 6: Time evolution of the density matrix element 
p e i (1,-1) for T = 100 mK and Q = 100. A) the modulus 
j/5 e i(l, —1)| and B) the phase normalized by n. 



cessional movement around the e y axis, that is, as L e i(t) 
oscillates between \L Z ) — | ± l) z (Figs. Q]and[H]) it passes 
through \L X ) = | ± l) x = (| + l) z ± | - 1),) /V2. This 
coherent dynamics persists during the relaxation time 
r e i w 100 ns and after that the rotational wavepacket 
looses most of the coherence. Furthermore, it evinces 
that the el-ph coupling with the odd parity phonon 
modes is the main responsible for the observed Rabi os- 
cillations. Finnaly, we point out that the orbital angu- 
lar momentum of the suspended QD system behaves like 
a spin 1/2 system and the interaction with the phonon 
modes of the platform play the role of a transverse mag- 
netic field. 



Electron-phonon dynamics in the presence of a weak 
magnetic field 

The dynamics of the coupled electron-phonon system 
is substantially changed in the presence of a magnetic 
field B applied perpendicularly to the plane of the QD. 
As noted in section , the diamagnetic energy term in Eq. 
(|4]) can be disregarded if the magnetic field is sufficiently 
weak. To satisfy the weak field condition, we set B = 

500 G, SO that Ediam/ 'E Zeeman = 0.1 and E d iar,i/E GS = 

0.008, where E GS is the energy of the electronic ground 
state | k) = |0, 1). In this regime the electron states with 
±Z are Zeeman splitted but the electronic wavefunctions 
remain unaltered. Thus, the ensuing effects have a purely 
kinetic origin, since they arise exclusively because of the 
elimination of the ±Z degeneracy. 

Figure [7] presents the dynamics of the electron-phonon 
system in the presence of B, disregarding mechanical 
losses for the moment, at the initial bath temperature 
T = 100 mK. For the sake of completeness, we also show 
the evolution of the initial state |^(0)) with L e /(0) = —1. 
The physics evinced in Fig. [7J should be contrasted with 
that presented in Fig. 0] 
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FIG. 7: Time evolution of the coupled electron-phonon state 
in the presence of a weak magnetic field, assuming only inter- 
nal decoherence and no external mechanical dissipation. From 
top to bottom the graphs correspond to: A) average electronic 
angular momentum L e i(t) = Tr e ;{Lp e ;(t)}, B) electronic pu- 
rity r e i(t) = Tr{/3gj(t)} and C) electronic energy. The curves 
correspond to initial states with L e i(t = 0) = ±1 and T = 
100 mK. The magnetic field is B = 500 G. 
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FIG. 8: Time evolution of the coupled electron-phonon state 
in the presence of a weak magnetic field, including attachment 
loss dissipation with Q = 100. A) average electronic angular 
momentum L e i(t) — Tr e i{Lp e i(t)} and B) electronic purity 
F e i(t) = Tr{pli(t)}. The curves correspond to initial states 
with L e i(t = 0) = ±1 and T = 100 mK. The magnetic field is 
B = 500 G. 



Two different physical effects occur for L e ;(0) = ±1, 
because of the Zeeman splitting of the initial state. In the 
case L e i(0) = — 1 the oscillations of L e i are almost com- 
pletely eliminated and T e i(t) w 1 for all times. The rea- 
son is that < lu±i < to a —x, rendering the state I = — 1 
decoupled from the rest of the system. As the magnetic 
field decreases, the oscillations of L e i quickly build up 
to the original value. The magnetic field also reduces 
the effect of the temperature in the el-ph relaxation, so 
that the electronic decoherence and the energy dissipa- 
tion processes are delayed, Figs. [7jB and C. 

For the case L e i(0) = +1, in addition to the decou- 
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pling < oj±i < w a= i, the Zeeman splitting produces 
the accidental resonance E(l, 1) — E(0,1) ~ huj a= n. 
That resonance sets up a strong Rabi oscillation between 
the I = +1 and the I — electronic states, which is 
responsible for the oscillation of L e i and, specially, the 
intense energy exchange with the phonon mode a = 17. 
The resonant effect is very sensitive to the detuning S. 
For the calculations presented in Fig. [7] we find that 
S = 0.004 x w Q= i7. However, by setting S — 0.02 x uo a= n 
the oscillation amplitude of L e ; decreases to a third. 

For a QD located at the center of the platform, the 
phonon modes of odd parity are the ones that generally 
couple more strongly with the electronic states. 

The mechanical dissipation affects the el-ph dynamics 
of the L e i — ±1 cases differently, as shown in Figure [8] 
Because the state with L e i(0) = — 1 exchanges very little 
energy with the phonon bath, the latter remains close 
to the thermal equilibrium and the attachment loss pro- 
duces a minor dissipation effect. The L e i(0) = +1 state, 
otherwise, exchanges a large amount of energy with the 
high energy phonon modes of the bath (Fig. [7]-C). Conse- 
quently, for those modes (n a ) 3> n^(T), which drives the 
phonon system out of equilibrium and leads to a strong 
mechanical energy dissipation. 



qubit structures, which consist of two adjacent quantum 
dots coupled by a tunneling mechanism, as well as from 
the spin-qubit devices. The proposed structure, never- 
theless, has characteristics in common to both original 
concepts, since the orbital angular momentum of the QD 
plays the role of the spin. 

Despite the interesting properties, the proposed NEMS 
poses significant challenges compared to the usual can- 
tilever and bridge structures. It is, therefore, desirable 
to consider alternative geometries for the platform, with 
the purpose of minimizing the attachment losses without 
removing the electronic oscillations. The current model 
evidences the prototypical characteristics of the struc- 
ture. 
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CONCLUSIONS 

We have investigated the quantum dynamics of the 
electron-phonon system for an unconventional design of 
nanoelectromechanical resonator. The interplay between 
well defined orbital angular momentum states of a circu- 
lar QD and the phonon modes of a wide nanomechan- 
ical resonator (here called platform) gives rise to a co- 
herent electron-phonon dynamics. The phonon ensemble 
is comprised by the quantized vibrational modes of the 
suspended platform and is described by a thermal field. 
Virtual phonon transitions cause the electronic state to 
undergo Rabi oscillations between the degenerate angular 
momentum states I = ±1, without direct energy transfer 
to the phonon ensemble. That fact circumvents the dis- 
sipation effects caused by the mechanical (attachment) 
losses in the nanoresonator. By analysing the electronic 
density matrix, it is evidenced that the orbital angular 
momentum of the QD behaves like a spin 1/2 system 
in the presence of a transverse magnetic field. It is also 
shown that a weak magnetic field can be used to suspend 
the dynamics of the electron-phonon system, by decou- 
pling the angular momentum states ±1 because of the 
Zeeman splitting. Several energy loss mechanisms are 
considered for the platform and the attachment loss is 
regarded as the main source of dissipation. 

The results cast light upon the underlying physics of a 
system that behaves like a single-QD charge-qubit, with 
the orbital angular momentum of the electron as the 
quantum bit variable. It differs from the usual charge- 



[5 

[7; 

[8 
[9 

[10; 

[11 
[12 

[13 
[14 

[is; 

[16 



P. Hawrylak, Phys. Rev. Lett. 71, 3347 (1993). 

C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 
L.P. Kouwenhoven, D.G. Austing, and S. Tarucha, Rep. 
Prog. Phys. 64, 701 (2001). 

J.M. Elzerman, R. Hanson, J.S. Greidanus, L.H.Willems 
vanBeveren, S. De Franceschi, L.M.K. Vandersypen, S. 
Tarucha, and L.P. Kouwenhoven, Phys. Rev. B 67, 
161308(R) (2003). 

A.K. Huttel, S. Ludwig, H. Lorenz, K. Eberl, and J. P. 
Kotthaus, Phys. Rev. B 72, 081310(R) (2005). 
J. Gorman, D.G. Hasko, and D.A. Williams, Phys. Rev. 
Lett. 95, 090502 (2005). 

T. Hayashi, T. Fujisawa, H.D. Cheong, Y.H. Jeong, and 
Y. Hirayama, Phys. Rev. Lett. 91, 226804 (2003). 

D. Loss and D.P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

J.R. Petta, A.C. Johnson, J.M. Taylor, E.A. Laird, A. 

Yacoby, M.D. Lukin, CM. Marcus, M.P. Hanson, and 

A.C. Gossard, Science 309, 2180 (2005). 

W.G. van der Wiel, M.Stopa, T. Kodera, T. Hatano, and 

S. Tarucha, New J. Phys. 8, 1 (2006). 

T. Brandes.Phys. Rep. 408, 315 (2005). 

A. N. Cleland, Foundations of Nanomechanics (Springer- 

Verlag, 2002). 

R. H. Blick, A. Erbe, L. Pescini, A. Kraus, D. V. Scheible, 
F. W. Beil, E. Hoehberger, A. Hoerner, J. Kirschbaum, 
and H. Lorenz, J. Phys.: Cond. Matt. 14, R905 (2002). 
A. D. Armour, M. P. Blencowe, and K. C. Schwab, Phys. 
Rev. Lett. 88, 148301 (2002). 

A. N. Cleland and M. R. Geller, Phys. Rev. Lett. 93, 
070501 (2004). 

L. G. C. Rego, A. Gusso, and M. G. E. da Luz, J. Phys. 
A 38, L639 (2005); A. Gusso, M.G.E. da Luz, and L.G.C. 
Rego, Phys. Rev. B 73, 035436 (2006). 



10 



X.M.H. Huang, C.A. Zorman, M. Mehregany, and M.L. [31 
Roukes, Nature 421 496 (2003). 

A. Gaidarzhy, G. Zolfagharkhani, R.L. Badzey, and P. [32 
Mohanty, Phys. Rev. Lett. 94, 030402 (2005). 
K.Y. Yasumura, T.D. Stowe, E.M. Chow, T.Pfafman, [33 
T.W. Kenny, B.C. Stipe, and D. Rugar, J. Microelec- 
tromechanical Systems 9, 117 (2000). [34 
P. Mohanty, D.A. Harrington, K.L. Ekinci, Y.T. Yang, [35 
M.J. Murphy, and M.L. Roukes, Phys. Rev. B 66, 085416 
(2002). 

X.M. Huang, X.L. Feng, C.A. Zorman, M. Mehregany, [36 
and M.L. Roukes, New J. Phys. 7, 247 (2005). 
M.O. Scully and M.S. Zubairy, in Quantum Optics (Cam- [37 
bridge University Press, 2002). 
PA. Cain, H. Ahmed, and D.A. Williams, J. Appl. Phys. [38 
92, 346 (2002). 

A.J. Fergusson, D.G. Hasko, H. Ahmed, and D.A. [39 
Williams, Appl. Phys. Lett. 82, 4492 (2003). 
R. Stomp et al, Phys. Rev. Lett. 94, 056802 (2005). [40 
P. Fallahi et al , Nanoletters 5, 223 (2005). 

G.J. Delben, A. Gusso, and M.G.E. da Luz, Phys. Rev. [41 
A 73, 055601 (2006). 

Karl A. Graff, Wave Motion in Elastic Solids (Dover, [42 
NY, 1975). 

D. Zhou, Y.K. Cheung, F.T.K. Au, and S.H. Lo, Int. J. [43 
Solids Struct. 39, 6339 (2002). 

A. Gusso and L.G.C. Rego, Phys. Rev. B 75, 045320 
(2007). 



H. Tal-Ezer and R. Kosloff, J. Chem. Phys. 81, 3967 
(1984). 

V.V. Dobrovitski and HA. De Raedt, Phys. Rev. E 67, 
056702 (2003). 

E.T. Jaynes and F.W. Cummings, Proc. IEEE 51, 89 
(1963). 

S. Sachdev, Phys. Rev. A 29, 2627 (1984). 
In the calculations we use C D p = 9.8 x 10 3 meV and 
m e = 0.95mo for the in-plane effective mass of the elec- 
tron. 

P.A. Greaney and J.C. Grossman, Phys. Rev. Lett. 98, 
125503 (2007). 

Y. Jimbo and K. Itao, J. Horolog. Inst. Japan 47, 1 
(1968). 

M.C. Cross and R. Lifshitz, Phev. Rev. B 64, 085324 
(2001). 

D.M. Photiadis and J.A. Judge, Appl. Phys. Lett. 85, 
482 (2004). 

M. Imboden, P. Mohanty, A. Gaidarzhy, J. Rankin, and 

B.W. Sheldon, Appl. Phys. Lett. 90, 173502 (2007). 

R. Lifshitz and M.L. Roukes, Phys. Rev. B 61, 5600 

(2000) . 

G. Zolfagharkhani, A. Gaidarzhy, S.B. Shim, R.L. 
Badzey, and P. Mohanty, Phys. Rev. B 72, 224101 (2005). 
D.T. Simon and M.R. Geller, Phys. Rev. B 64, 115412 

(2001) . 



